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Abstract 

This report summarizes centrifuge-related work performed at the Smart Systems Research Laboratory at NASA 
Ames Research Center’s Computational Sciences Division from 1995 through 2003 The goal is to develop an 
automated system that will sfense an imbalance (both static and dynamic 3 ) in a centrifuge and issue control commands 
to drive counterweights to eliminate the effects of the imbalance. This autobalancing development began when the 
ISS centrifuge design was not yet finalized, and was designed to work with the SSRL Centrifuge laboratory 
prototype, constructed in 1993-1995 and shown in Figure 1 Significant differences between that prototype and the 
current International Space Station (ISS) Centrifuge design are that 

• The spin axis for the SSRL Centrifuge prototype can translate freely in x and y, but not wobble, whereas the 
ISS centrifuge spin axis has 3 translational and two rotational degrees of freedom, supported by a vibration 
isolation mechanism. 

• The imbalance sensors are strain gauges both in the rotor and the stator, measuring the imbalance forces, 
whereas the ISS centrifuge uses eddy current displacement sensors to measure the displacements resulting 
from imbalance. 

High fidelity autobalancing and FDIR systems (for both counterweights and strain gauges) are developed and tested 
in MATLAB simulation, for the SSRL Cmtrifuge configuration. Hardware implementation of the autobalancing 
techno logyjfeegun in 1996, but was terminated due to lack of funding. The project lay dormant until 2001-2002 when 
the FDIR capability was added, 

Now that the IS S centrifuge configuration is (close to) finalized, the logical extension of this research is to extend its 
■ application to the new configuration to evaluate the potential autobalancing performance improvements and FDIR 
capabilities, Following are the steps towards completing this task, as of 25 July 2003 

1 Removing the constraint that the spin axis is fixed, allowing free or passively/actively constrained motion of 
the spin-axis. 

2. Represent imbalance differently, to account for 3-D model - e.g., mass center location and inertia matrix 

3 . Update to use the present counterweight locations. 

4 Update to use the present bearing displacement sensors (BBSs) - this would seem to force me to undo the 
step where all sensors are first reduced to a net force and torque on the rotor 


1 Ed.Wilson@intelUzation.com, 454 Barkentine Lane, Redwood Shores, CA 94065-1126 

2 RobertW.Mah@nasa.gov, SSRL group lead 

3 Here, in accordance with standard rotor-balancing terminology, a “static” imbalance is one that creates a net force 
•or torque on the centrifuge rotor in the presence of gravity when the rotor is not spinning (i.e. the center of mass is 
not on the axis of rotation). This is the reason for the name “static”, although there will be no imbalance force created 
in zero g with the rotor stopped. A “dynamic” imbalance creates a net force or torque on the centrifuge rotor when 
the rotor has a non-zero angular velocity or angular acceleration (i.e. ; tire cross-axis terms in the inertia matrix are 
non-zero). 
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5. Model the vibration isolation mechanism (VIM), including springs and active dampers 

6. Integrate with ADAMS for testing - thought is to have a simpler possibly linearized model derived (by 
hand) for the autobalance^ design, but then to test it on the ADAMS model that is fully nonlinear and models 
the rotor flexibility 

Table of Contents 

1 Introduction. 6 

11 SSRL Centrifuge laboratory prototype. 6 

1 .2 Brief summary of the autobalancing method presented in this report: 1 

2 Control System Architecture, U 

3 Imbalance Identification 11 

3 . 1 Description of a structure to model a general state of imbalance in the rotor (defined by a set of “model- 

imbalance parameters”). 12 

3.2 Derivation of equations relating model-imbalance parameters to the net forces and torques on the rotor 
16 

3.2 1 Additional forces 20 

3.2.2 Summary 23 

3 .3 Derivation of equations relating counterweight parameters to the net forces and torques on the rotor 24 

3.4 Identification of the net forces and torques based upon the strain gauge signals. 28 

3 . 5 Calculation of forces due to imbalances alone by subtracting out the effects of the counterweights. 37 

3.6 Identification of the model-imbalance parameters based upon these net forces and torques. 39 

3 7 SG fault detection and identification 44 

3.8 A MATLAB simulation to validate the identification scheme, including sensitivity to noise in the strain 

gauge readings 45 

3.9 Filters to produce \f : ffi, and a from the tachometer and encoder signals. 49 

3 .10 Filters to produce F( A +b) x , F(A+B) y > F(a-b>, and F( A -B) y from the strain gauge signals . 50 

4 Calculation of the desired counterweight positions based upon the identified imbalance 53 

5 Servo control loop to move the counterweights to their desired locations. 54 

6 Development approach ' Work plan. 54 

Table of Figures 

Figure 1 The S SRI Centrifuge Autobalancing Hardware S imulator . 6 

Figure 2 Control system block diagram U 

Figure 3 Coordinate system used to describe imbalances and their effects 1 3 

Figure 4' z-axis locations of strain gauges, counterweights, and imbalances 14 

Figure 5 : Model imbalances . 1 5 

Figure 6: Misalignment between the axis of rotation and the gravity vector results in additional forces that are not due 
to imbalances. 21 

Figure 7- Counterweight coordinates. .. .. 25 

Figure 8: Strain gauge layout. 29 

Figure 9: Software simulation GUI - bearing forces 46 

Figure 10; Software simulation GUI - imbalance tracking, Strain gauge FDIR 47 

Figure 1 1 . Software simulation GUI - simulation control panel 48 

Figure 12 Software simulation GUI - strain gauge outputs. . 48 


SRPT f'e.ntrifUaf*. AivlohaVanr.tricr and Pr>TP 


PdwarH Wilson 


75 Tnlv 7nni 


naocp. 7/70 



1 Introduction 


1.1 SSRL Centrifuge laboratory prototype 

The SSRL Centrifuge laboratory prototype is shown in Figure 1 It was developed in 1993-1995 by Robert Mah and 
Mike Guerrero at the SSRL, with Alessandro Galvagm assisting with the computer systems interfaces. It floats on 
four 8-inch diameter Fox Air Bearings with spherical air bearing pivots, allowing the vertical spin axis to translate in x 
and y in response to imbalances. There is slight axial motion permitted (see below), but it is a much stiffer suspension 
than the VIM of the ISS Centrifuge. 



Figure h The SSRL Centrifuge Autobalanting Hardware Simulator 


There are 12 counterweights (CWs) constructed of stepper motors (providing the counterbalance mass as well as 
actuation) and encoders that drive along lead screws (also known as ACME screws or trapezoidal screws these are 
different from the ISS Centrifuge which uses ball screws). The CWs are located in upper and lower planes of 4 radial 
CWs and 4 vertically moving CWs. The vertically moving ones are not used in this research, and the 8 remaining 
CWs provide redundancy as only 4 are required. 

There are locations for 24 strain gauges (SGs) to measure imbalance forces, 8 do not rotate, and are for measuring 
the force between the stator and the air bearing platform. No SGs are installed at these locations. 8 SGs rotate with 
the rotor, near the upper CW plane. Of these 8 locations, 4 SGs are installed, equally spaced through one rotation. 
Same goes for a plane near the lower CW plane, meaning that the hardware presently has a total of 8 rotating SGs 
The SG assemblies are visco-elastic, allowing for slight displacement (on the order of mm) and providing damping. 
The angular (tip-tilt) natural frequency of the rotor is on the order of 3 Hz, and the damping ratio is about 0.1 
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Planned, but not yet installed, are displacement sensors that would measure the base translations in the horizontal 
dimensions. 

1.2 Brief summary of the aiitobalancing method presented in this report; 

1 Centrifuge rotor is modeled as a rigid body spinning about a fixed axis. 

2 Imbalance is represented in a compact, intuitive way as the x- and y-iocations of a point mass in the central 
plane and a pair of asymmetrically located point masses in off-central planes. This four-parameter 
representation is sufficient to represent an arbitrary imbalance and can be intuitively related to counterweight 
motions. 

3 The rotor imbalance (not including counterweights) is estimated at each sample period as follows: 

a. Sensor signals (strain gauges, counterweight positions, velocities, and accelerations, rotor angle 
encoder and tachometer) are combined, using a least-squares fit, to calculate the estimated net 
force and torque in x and y (F^, FyiB, ^xib> t y m) created by the imbalance 

b. The four imbalance parameters (mentioned in (2), above) and their derivatives are estimated using 
these forces and torques This uses a dynamic model of the rotor that calculates, effects due to the 
position, velocity, and acceleration of the imbalance parameters 

4. The counterweights are driven to exactly counteract the estimated rotor imbalance. 


This can be considered an indirect method (analogous to indirect vs, direct adaptive control) since the sensor signals 
are used to build a model of the imbalance, then corrective action is taken based on the identified model parameters 
In a direct method, (filtered, and mathematically manipulated) sensor signals would be used to directly drive the 
counterweights. Hopefully, this feedback loop would drive the counterweights until the sensors read zero. The 
increased complexity of the indirect method presented here enables more accurate fitting of the sensor data to the 
dynamic model of the imbalance. Whether the increased accuracy of the indirect method produces results that are 
sufficiently better than those of the direct method will depend on the characteristics of the imbalance (how fast it is 
moving, etc.) and sensor noise. 

A prior approach was developed by the author, concluding in October 1995, that was very similar The major 
improvement made in this updated version involves breaking up the identification so that all sensor signals are 
reduced (combined) to result in an intermediate estimation of the net forces and torques on the rotor In the previous 
approach, all sensors were used directly to identify the imbalance parameters. In the present approach, all sensors 
(strain gauges, counterweight positions, velocities, and accelerations, rotor angle encoder and tachometer) are used at 
each sample period to calculate four variables; the estimated net force and torque in x and y (F^, F>ib, ^ib) 
created by the imbalance. This process involves a least squares fit to the data, using a model of the sensor geometry 
and subtracting out known forces due to the counterweights. These four variables then pass to the imbalance 
identificationalgorithm. Benefits of segmenting the identification into these two parts are: 

1 Physically, the four imbalance parameters are directly related to these four. intermediate variables (Fxib, 
t*©, The relation between sensor values and these intermediate variables is more direct than that 
between sensor values and imbalance parameters. This logically separates estimation of forces and torques 
created by the imbalance from the estimation of imbalance parameters themselves. 

2 It is easier to identify failed sensors, since the analysis can be performed without regard to the imbalance 
dynamics - one can analyze the residuals in the estimation of the intermediate variables. 

3 . If sensors change (e.g.. on-line failure, design change, etc.), the second part of the identification (that finds 

imbalance parameters from estimated forces and torques) does not have to be changed 

4. Overall complexity is reduced by breaking one large problem into two smaller ones. No accuracy is lost, due 
to the physical reasoning listed in (1) above. 
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2 Control System Architecture 



ID system: use sensed variables 

(T.® & V -V 

and known parameters to estimate 
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Figure 2: Control system block diagram 


3 Imbalance Identification 

The basic approach taken here to identify the imbalance is: 

1) Calculate (identify) the net forces and torques on the rotor at each sample period, based upon the strain gauge 
measurements. 

■ - 2) Subtract out the forces due to the counterweights, leaving the forces due only to the imbalances . 

3) Identify the “imbalance parameters” corresponding to these forces that define the state of imbalance in the rotor. 

An indirect, model-based approach like this should work well if the form of the model can be identified correctly and 
the sensors are not excessively noisy or biased. 

3.1 Description of a structure to model a general state of imbalance in the rotor (defined by a set of 
“ model-imbalance parameters”). 

The first step here is to select a general form for the imbalance that is capable of representing any possible imbalance. 
The 3x3 inertia matrix describing the spun portion of the centrifuge (the “rotor 15 ) could be used, but a more intuitive 
representation is presented here. This arbitrarily chosen imbalance model structure has been proven to be sufficiently 
general to model any possible imbalance. 4 

The form chosen contains two specific perturbations to a perfectly balanced rotor These “model imbalances” are: 

1 ) A “point-mass” imbalance (PMI) located on the plane equidistant from the two counterweight planes. 


4 See Appendix C 
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2) A “mass-couple” imbalance (MCI) composed of two equal masses located symmetrically about the center of the 
coordinate system. Each mass is in one of the planes containing the counterweights. 

It is not possible to represent a general state of imbalance with only a single point mass located somewhere within the 
rotor This is most easily demonstrated hy considering a perfect rotor that has two equal masses added to it 
symmetrically about the mass center (similar to the “mass-couple” imbalance above). This situation can not be 
represented by a single point mass imbalance, proving that it is not a general representation. 


The forces acting on each of these model imbalances (PMI, MCI) are due to gravity, rotor motion (centrifugal 
acceleration and angular acceleration), and imbalance motion (Coriolis effects, imbalance acceleration) 3 * 5 The net 


forces and torques [F x , F y , T*, T y ] on the rotor will be calculated as functions of the rotor motion (\|/ 3 oo, a), fixed 
system parameters (g, rotor geometry), and imbalance parameters. Once measurements of actual net forces torques 
and rotor motions are available, it will be possible to estimate the imbalance parameters. 



x, y, zare fixed In the disk frame 
X, Y, 2 are fixed on the de-spun portion of the 
centrifuge 

is the angle of the disk the only rotationa 
degree of freedom 

A, B represent the upper and lower portions of 
the disk 

z = 0 at the midpoint between the 
counterweight planes 


Figure 3: Coordinate system used to describe imbalances and their effects 

The coordinate system is shown in Figure 3 ; with the z-axis locations shown in detail in Figure 4 An xyz Cartesian 
coordinate system is used to describe imbalance locations fixed within the rotor The rotor angle, \|/( t), describes the 
relative angular position of the rotor with respect to the non-spinning (i.e. U de-spun”) portion of the centrifuge 6 The 
rotor angular velocity, co(t), and angular acceleration, a(t), are the first and second derivatives of V|/(t) with respect to 
time. A complete list of variables is contained in Appendix A. 


3 The rotor angular acceleration force, Coriolis force, and force due to imbalance acceleration are expected to be 

minimal, but will be included for completeness. 

6 Although the centrifuge base may actually rotate slightly, this will be a small effect, and it is neglected in this 

analysis. 
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A (B) refers to the upper (lower) part of the disk 
m p , m c are the imbalance masses 
m is the mass of the disk 

CW A . CW B are the planes of the counterweights (z = +/- 1 by 
definition) 

S As is the plane of the upper (A) spinning (s) strain gauges (S) 
S Af is the plane of the upper (A) fixed (f) strain gauges (S) 

$ Bs is the plane of the lower (B) spinning (s) strain gauges (S) 
$ sf is the plane of the lower (B) fixed (f) strain gauges (S) 
z m is the z-axis coordinate of the disk c.o.g. 

is the z-coordinate of S As 
z Af is the z-coordinate of S Af 


is the z-coordinate of S R 


z Bf is the z-coordinate of S Bf 

z - 0 is chosen arbitrarily as the midpoint between CW A and CW B 


Figure 4: z-axis locations of strain gauges, counterweights, and imbalances 


The origin of the xyz coordinate system is fixed on the axis of rotation at the midpoint between the upper (A) and 
lower (B) counterweight planes. The “mass-couple” imbalances, shown in Figure 4, are defined arbitrarily to lie in the 
same planes as the counterweights. The ‘-point-mass” imbalance is chosen arbitrarily to lie in the plane equidistant 
between these two planes. The z-axis coordinate of this plane is defined to be z = 0 The z-axis coordinates of the 
upper and lower counterweight planes are defined to be z = f and z - -f 

Two types of strain gauges will be used on both the upper and lower parts of the rotor “Fixed” gauges (Sat and SBf) 
will be attached to the non-spinning base, and “spinning” gauges (Sas and S Bs ) will be mounted on the spinning rotor 
As their z-axis coordinates may not be symmetric, each of the four strain gauge locations has its own coordinate 
label, Zas, ZAf, z B f> and z Bs , as shown in Figure 4. 

To summarize, the physical locations of the counterweight planes define the locations of the z-axis and the model 
imbalances (PMI and MCI). 

F a is the force on the point masses 
due to angular acceleration (not 
shown in the figures) 

F c is the force on the point masses 
due to centrifugal acceleration 
F g is the force on the point masses 
due to gravity 

F x is the net force on the disk in the 
+x direction (F y is not shown) 

X y is the net torque about the y axis 

of the disk (T x is not shown) 

Forces due to imbalance motion 
(cpriolis and acceleration) are 
not shown here. 




Figure 5: Model imbalances 

Point-mass imbalance (PM}> Mass-couple imbalance (MCI), and their effects on the rotor 


Figure 5 shows the model imbalances used to describe the general state of imbalance in the rotor These discrete point 
masses, if placed in the correct locations, can produce the same effects as any imbalance in the rotor This is proven in 
Appendix C and summarized as follows. The proof assumes that any state of imbalance can be described by a 
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(possibly large) number of point masses of varying mass and location. The net forces and torques due to each of these 
many masses is calculated, summed, and shown to be equal to the net forces and torques due to the two model 
imbalances (PMI, MCI), provided the imbalance parameters are chosen correctly (for example, m p x p = S i=!l n x w 
where mi and x* represent the mass and x-coordinate of each of the many point masses). In Section 3.2, the effective 
net forces and torques on the rotor due to the model imbalances will be calculated, 

3,2 Derivation of equations relating model-imbalance parameters to the net forces and torques on the 
rotor, 

Even though the spin axis is not perfectly rigid, it is modeled as such in this analysis. 

The first model imbalance, the point-mass imbalance, has a mass equal to m p and is located in the x-y plane, with z = 

0 The mass magnitude, x-y location, x-y velocity, and x-y-z acceleration all affect the net forces and torques on the 
rotor in the x and y directions (the effect due to acceleration in the z-direction is small, and will later be dropped). 
Therefore, [m p , x p , y p , v xp , v >p , a xp , ayp] define the imbalance parameters for the PMI (the value, of m p is arbitrary). 

Note that these forces and torques are computed in the rotor frame. For example, F x is the net force on the rotor 
aligned with the -f-x axis of the rotor. Since forces and torques aligned with the z-axis will not be measured by the 
strain gauges, they are not included here, The point-mass imbalance results in. the following net forces and torques: 

• Fg, the force due to gravity. This acts , straight down at all times 7 and is independent of rotor motion. Its effect on 
the rotor appears as a torque equal to m p gr p , where r p is the radial distance from the axis of rotation. 

= 0 F y = 0. X x ~ -m p gy p T y = m p gx p 

• F c , the centrifugal force. This acts radially outwards whenever the rotor is rotating, and is proportional to the 
angular velocity squared, Its effect on the rotor appears as a force in the direction of the point mass location 
vector (x Pj y p ) and with a magnitude equal to m p r p CQ 2 

F x ~ m p co 2 x p F y = m p co 2 y p T x = 0 T y ~Q 

• F a , the force due to rotor angular acceleration. This force acts perpendicular to the point mass location vector 
(x p> y p ) whenever the rotor has an angular acceleration, and is proportional to the angular acceleration. Its effect 
on the rotor is a force in the direction perpendicular to the point mass location vector and with a magnitude equal 
to m p r p a. 

F x “m p ay p F y ~ -m p ax p T x ~0 T y = 0 

• F v , the Coriolis force due to the imbalance moving within the rotor This force acts perpendicular to the 
imbalance velocity, (v xp , Vyp) and is proportional to the imbalance velocity and the rotor angular velocity. 

F x “2m p ©v yp F y = -2m p ©v X p T x - 0 T y ~0 

• F a , the force due to the imbalance accelerating within the rotor This acts in toe opposite direction of the 
acceleration, (a xp , a^, a^). 

F x — un p a xp Fy m p ayp T x 1 "m p a 2p yp T y m p a.£pX p 


The second model imbalance, the mass-couple imbalance, consists of two point masses of mass m c . The z-axis 
location of the upper (lower) point mass is chosen arbitrarily to lie in the same plane as the upper (lower) 
counterweights. Since the masses are located symmetrically about the geometric center of the rotor by definition, the 
x-y location of the upper point mass (x c , y c ) also determines the lower point mass location (-x c> -y c ). There is no net 
force due to gravity, also due to symmetry Positions (x c> y c ), velocities (v^ v JC ), and accelerations (a xc , a^) of the 


7 Effects due to misalignment between the gravity vector and the axis of rotation are described in Section 3 .2. 1 
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imbalance masses are always equal and opposite by definition. Therefore, [m c> x c , y c , v xc , v JtJ , a xc , a } J define the 
imbalance parameters for the MCI The MCI results in the following forces on the rotor in the rotor frame, 

• F g , the force due to gravity, is zero by definition. This is because the masses are defined to be symmetrically 
located about the center of the rotor coordinate system. 

F*=0 Fy = 0 T X ™ 0 Ty=0 

• F Cj the centrifugal force. This acts radially outwards whenever the rotor is rotating, and is proportion^ to the 
angular velocity squared. 

F x - 0 F y ~ 0 X x = -2mc(fi 2 /y c T y - 

• F a , the force due to rotor angular acceleration. This force acts perpendicular to each point mass location vector 
whenever the rotor is accelerating, and is proportion^ to the angular acceleration. 

F x = 0 F y = 0 T x = 2m c a/x c T y =2m c a/y c 

• F v , the Coriolis force due to the imbalance moving within the rotor: This force acts perpendicular to the 
imbalance velocity, (v** v )c ) and is proportional to the imbalance velocity and the rotor angular velocity. 

F x — 0 Fy~0 T x = 4m c ©/v xc T y = 4ro c co/v >c 

• F fl , the force due to the imbalance accelerating within the rotor This acts in the opposite direction of the 
acceleration, (a xc , a yB> a zc ). 

F x = 0 F y ~0 X x = 2m c /a jc 2m c a zc y c T y =*= -2m/a xc + 2m c a zc x c 


Adding the net forces and torques on the rotor that are created by the two imbalances (forces are calculated in the 
rotor frame); 



m p © 2 Xp 

+ m p ay p 

+ 2m p a>Yyp 

m p a xp 



Fy= 

m p © 2 y p 

m p ax p 

2m p ©v xp 

m p ayp 



Tx = 

• m pgy P 

mp^zpYp 

2m c a> 2 /y c 

+ 2m c a fx<. 

+ 4m c a>/v xc + 2mc/a yx . 

imc^yc 

T y = 

m p gx p 

+ m p a zp x p 

+ 2m c a> 2 /x c 

+ 2mca/y c 

+ 4m c ©/v JC 2mA® 

+ 2m c a.,, c x c 

Putting the equations in a linear form so that the linear regression can be applied later requires dropping the “cross- 
product” terms. These are some of the terms due to acceleration of the imbalance within the rotor frame likely to be 
small. # 


m p C0 2 Xp 

+ m p ay p 

+ 2m p coVj P 




Fy = 

m p ca 2 y p 

m p ax p 

2m p G)V xp 

m p a w 



t:*= 

-m p gy p 

2m c co 2 /y c 

+ 2m c a/x c 

+ 4mc m/v xc 

2mc/a yp 


ty = 

m p gx p 

+ 2m c CD 2 /Xc 

+ 2mca/y c 

+ 4mcC o/vjp 

2m c /a xc 



3.2.1 Additional forces 

The misalignment between the rotor axis of rotation and the gravity vector will cause strain gauge measurements that 
are not due to the imbalances. These forces must be accounted for so that the model imbalance parameters may be 
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properly identified. The misalignment is quantified by the misalignment angle in the x and y -directions, <f> x and <j> y , This 
misalignment causes a constant force on the fixed strain gauges and an oscillating forces on the spinning strain 
gauges. 

The misalignment causes forces due to the model imbalances as well, but as these are second order effects, they will 
not be considered. These terms will not apply for the space station centrifuge, as it is in a micro-gravity environment 



Figure 6: Misalignment between the axis of rotation and the gravity vector results in additional forces that 

are not due to imbalances 


F s ^, the force due to the axis misalignment, is affected by the total mass of the rotor, m, the z-axis location of the 
rotor center of gravity (c.o.g.), Zj„, and the misalignment angles, (j> x and Angles are assumed to be small, so sin(<|>) 
* <f>. It results in the following net forces and torques on the rotor measured in the base frame. 

Fx = mg<{> x Fy - mg<j>y T x ^-mg(() y z m X Y == mg<|> x z m 

Unlike the forces due to imbalances, these forces do not rotate with the rotor To calculate the forces in the rotor 
frame, they are rotated by an angle -ij r 

F x = FxPOs\[ / + F Y sin\}/ F y - -Fxsin\|/ + F Y cos\|/ T x = T x cos\|/ + tysinij/ X y = ■Txsinij/ + Tycosij / 

Substituting. 

F x = mg<j> x cosij/ + mg<fiy simp F y = *mg<j> x sinij/ + mg<}) y cosij/ 

T x = -mg<{) y z m cosvj/ + mg(j) x z nl sinij; T y = mg(|) y z m sinij/ + mg$ x z m cosy 

The parameters describing the imbalance must appear linearly in the equation if they are to be identified using a linear 
regression, so the torque terms containing z m <j> x and z m <j> y must be changed. New variables z m $ x ( = z m <|) x ) and z m $ y ( 
= z, n <j> y ) are introduced to maintain linearity. 

F x = mg(j> x cosij; + mg<j) y sin\)/ F y = -mg(j) x sinij/ + mg(|) y cosij; 

X x = -mg z nl4y cosij/ + mg z m+x sinij/ X y = mg z l4y sinij/ + mg z m!j>x cosij; 

Adding these forces and torques to the model- imbalance forces and torques, 

F x = m p co 2 Xp + m p ay p + ImpCov^ m p a xp + mg cosij/ (j> x + mg sinij/ <j> y 
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Fy = 

m p co 2 y p 

m p ax p 

2m p ©v xp 

rripSyp 

mg siny <}> x •+ mg cosy <j> y 


-m P gy P 

2m c <D 2 /y c 

+ 2mca/Xc 

+ 4m cCa/v xc 

+ Imc&yc + mg siny z^ x mg cos y z n ^ y 

^y = 

nipgXp 

+ 2m 0 co 2 /x c 

+ 2m c a/y c 

+ 4mcCo/v^ 

2mc/a xc + mg cosy z m ^ x + mg simp z m $ y 


Putting these equations in matrix form, 

FnetIB = 08 

where 0 is a 16x1 column vector of the model-imbalance parameters, m, m p , and m c may be chosen arbitrarily, since 
they simply scale the values of the other parameters. They are chosen to equal 1 (kilogram) so they drop out of the 
equations completely. is a 4x16 matrix, (sy s siny, cy * cosy). 


0 * 
x p 

Yp 

*c 

Yc 

Vxp 

v )p 

Vxo 

Yw 

a xp 

a JP 

3-XC 

a yc 



Zmijjx 

Z ™fo' 

<E> = 

CD 2 

a 

0 

0 

0 

2co 

0 

0 

1 

0 

0 

0 

gey 

gsy 

0 

0 

-a 

co 2 

0 

0 

-2co 

0 

0 

0 

0 

■1 

G 

0 

•gsy 

gey 

0 

0 

0 

*g 

2a f 

-2e?f 0 

0 

4 of 

0 

G 

0 

0 

2 f 

0 

0 

gsy 

-gey 

g 

0 

2 cz 2 f 

2c if 

0 

0 

0 

4o/ 

0 

0 

■if 

0 

0 

0 

gey 

gsy 


3.2.2 Summary 

The above linear matrix equation describes the net forces and torques on, the rotor, F net iB ~ Pv F y , T*, T y ] T that are 
due to the model-imbalance parameters, 0 = [x^ y p , x* y c , y xp> Vyp, v xc> v JC , a xp , a^, a** a^, ^ <j> y , z n4x , z m * y ] T which 
represent the model imbalances as well as the axis misalignment <j> y , z m ^ z m $ y ]. Note that the MCI produces 
torques only no forces, and the PMI produces no torques except for the gravity effects. 

Section 3.4 presents a linear regression method to estimate the model-imbalance parameters based upon force and 

torque measurements. However, sensors to measure [F x , F y , T x , T y ] directly are not available, so before this method 
can be applied, these signals must first be estimated based upon the strain gauge measurements. 

3.3 Derivation of equations relating counterweight parameters to the net forces and torques on the rotor. 

If the forces and torques due to the counterweights can be calculated and subtracted from the measured forces and 
torques, the remaining signal can be used to identify the imbalance alone. The basic approach to deriving the effects 
of the counterweights is to calculate the equivalent set of mo del- imbalance parameters and then use the equations 
already derived in Section 2b. 


Problem definition. 

Counterweight parameters, m p cw> m c cw, and 8ew= [Xpcw, y P cw, x c cw» Yccw> v xp cw, Vjpcw, Vxccw, v }t; cw, a xp cw, a ypC w, 
axccw, a >c cwl T define the counterweight properties that result in net forces and torques on the rotor These 
parameters are defined with the same structure as used for the model imbalances, 

Take measurements of counterweight positions, velocities, and accelerations, and calculate the above parameters. 
Then use these parameters to calculate the net forces and torques applied to the rotor by the counterweights 

This is to be performed at each sample 
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5 1t 8 2 coordinates describing position 
of upper moving counterweights 

8 3 , are in lower plane, not shown. 

x pcw- etc - are based upon 5 3 = 

8 n and 8 4 = 0. 

x pcw» Ypcw effect ' ve counterweight point 
mass location ( = (5 v B 2 /2) for 
the coordinates given) 

XcCW* y cCW effective counterweight mass 
\ couple location { = (0, 6 2 /2) for 

i the coordinates given) 


Figure 7 ; Counterweight coordinates 

The counterweights are arranged in the layout shown in Figure 7 The positions of each of the four counterweights is 
described by [Si, 8 2 , 6 3 , S 4 ]. 6j is the counterweight in the upper plane in the +x direction, 8 2 is the counterweight in 
the upper plane in the +y direction, 5 3 is the counterweight in the lower plane in the +x direction, S 4 is the 
counterweight in the lower plane, in 1 the +y direction. m C w is the mass of the moving portion of the counterweight 
assembly. The origins of the- axes for 8), S 2 , 5 3 , and S 4 are chosen so that with 5 = 0 , the net effect of the 
counterweight assembly is zero, regardless of where the counterweight actually is this is the position where the 
counterweight is exactly counterbalanced by the dead weight 

The velocities and accelerations are described by the first and second derivatives of the position signals, [ 5 j 8 2 .. 5 3 
5 4 §i S 2 8 3 ’ 84”]. Sensors will not be used to measure these signals directly (i.e., no tachometer feedback), and 
since the velocity and acceleration effects are likely to be small, they will not be calculated in the autobalancing 
system. However, their effects are included here for completeness. For this analysis, the counterweight coordinate 
vector, [8 l,8 2 , 8 3 , 8 4 , 8.i $ 2 S 4 §1” 8 2 S3" & 4 ”], is assumed to be available at all times. 

The first step is to calculate the counterweight parameters 0 C w - [Xpcw> y P cw> *ccw> y c ew> v xpC w> v Jp cw> v^ C w, v >cC w, 
axpcw: %ew> a x cC w, Vcw] T from the counterweight coordinates, [ 5 Ia 5 2 , 83, 5 4> 8! 5 2 5 3 8 4 Sj" S 2 ” S3” 84"]. 
The basic approach is similar to the proof shown in Appendix C add the individual effects of each of the 
counterweights and select the counterweight parameters to yield the same effects The equations are simplified if the 
counterweights all have mass, equal to mew (this- is the moving portion of the counterweights). 


%CW~ 

4 mew 


XpCW = 

(8-i + 83) / 4 


ypew^ 

(8 2 + 84) / 4 


VxpCW ” 

d/dt(x pCW ) 

— ( 5 j + 83 ) / 4 

VjpCW = 

d/dt(y p ew) 

~ (8 2 + 84 ) / 4 
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(6|" + 8 3 ")/4 
(Sj + 84 ) / 4 


a«pcw “ d/dt(v^ C w) 

aypcw = d/dt(v >pCW ) 

mccw = 2 m C w 

^ gcw ^ (81 -83) / 4 
y c cw ^ (82 - 84) / 4 

Vxccw= d/dt(x cC w) - (&! - 83 ) / 4 

v^cw = d/dt(y cC w) ~ (82 -- 5*)/ 4 

axccw = d/d^Yxccw) “ (81 - 83 >/ 4 

a>ocw= d/dt(Vj C cw) = (82 - 84 ) 7 4 

Substitute these counterweight parameters into the model-imbalance equations: 


Fxcw = 

2 

m p cwt° x pcw 

+ m p cw^y P cw 

+ 2m p cw®v yp cw 

mpcwaxpcw 

F>cw- 

-m pC w^ x pCW 

+ m p cw^ 2 y P cw 

2m pC w<av X p CW 

IppCWEljpCW 

TxCW = 

m P cwgy P cw 

+ 2mccwCt/Xccw 

2mccw© 2 /yccw 

+ 4m c cw®/VxcCw + 2m cC w^cCW 

T>CW ~ 

mpewgXpew 

+ 2mccwa> 2 /xccw 

+ 2mccw0t/y c cw 

+ 4m c cw®/Vj«cw 2m c cw/axccw 


Propping terms including counterweight velocity and acceleration 
Fxcw ~ ^pcw0> 2 ^pcw + m p cw«y P cw 

Fjcw = -m p cw ax pcw + m p cwC0 2 y p cw 

T x cw - m p cwgy P cw + 2mccwa/x cC w 2mccw© 2 ^cCw 

T?>cw ~ u^pcwg^pcw + Imccw^^ccw + 2m c p^afy c cw 

Putting these equations in matrix form, 

Fneicw — ^CW @CW 

where Qcw is a 4 x 1 column vector of the model-imbalance parameters O cw is a 4 x 4 matrix 
Qcw - [xpcw y P cw Xccw y c cw] T 

0 C w 855 nipcw® 2 m p cw^ 0 0. 

m p cw a m p cw <n 2 0 D 

0 m p cw & 2 mccw&f -2 m c cwt° 2 ^ 

m pC w g 0 2 mcew^r 2m c cw°d r 

3.4 Identification of the net forces and torques based upon the strain gauge signals. 

The goal in this section is to develop a procedure to produce net force and torque measurements, [F* F y , T x , T y ], at 
each sample period based upon strain gauge measurements. There will be between four and eight (and possibly more) 
strain gauge measurements used to produce the four net force and torque signals, Strain gauges maybe spinning with 
the rotor or fixed in the base Each strain gauge will have a bias, which must be accounted for The approach taken 
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here is to solve a least-squares fit to get the net forces and torques, while keeping the bias terms separate so that they 
can be identified in the imbalance identification step. 

For now, assume the following strain gauge layout: 

16 strain gauges are used 8 fixed, 8 spinning. 

They are arranged in sets of four. The upper set of fixed strain gauges is shown in the following sketch “S” indicates 
“strain gauge.” “A” or “B” (not shown here) indicates the upper or lower set. “x” or “y” indicate the axis of 
measurement w f * Or “s” (not shown here) indicate “fixed” or “spinning” gauges. “1” pr “2” identify gauges on the 
same axis, “1” is on the positive side of the axis and “2” is on the negative side of the axis. When there is a force on 
the rotor in the positive direction, “1” will be in tension and “2” will be in compression. 



Figure 8: Strain gauge layout 

The first step is to combine two strain gauges on the same axis, which should yield exactly opposite readings. To 
improve accuracy, these readings will combined by subtraction For example, = (Saxti Saxg)/ 2 If a force of -HO 
Newtons is applied to the rotor in the +x direction, Saxu will be in tension, reading -HON and Saxe will be in 
compression, reading TON The combined value, Sax& will be (IQ (• i0))/2= ION. This combined value will also be 
filtered to reduce sensor noise and extraneous vibrations For now, this operation is represented as “filterQ” 


^Axf ** filter! (Saxa Saxe) n) 

SAyf * filter( (S.Ayfi Save) / 2 ) 

Ssxf ** filter! (S Bx n S Bx e) / 2 ') 

Sfiyf " fxlter( (S B> fi Sb v e) / 2 ) 

Saxs “ filter! (Saksi Saxs 2 ) / 2 ) 

S Ays = filter! (SAysl SAysl) 1 % ) 

S Bx s - filter( (Sbxsi Szxst) 1 2) 

Sb^s ■=*■ filter( (S B )si 2 )/ 2) 


Assume, that this force measurement is composed of three parts: 

S - S.tn» + P + £ 

1) Stme, the actual force transmitted by the combined pair of strain gauges (if the sensor were perfect, S = S tme ). 
However, S lTue may contain force disturbances that are not due to imbalances (such as structural vibrations). 

2) p, a bias term that represents the net bias of the combined pair of strain gauges (assumed to be non-time- 
varying). With eight strain gauge pairs, there will be eight biases, 

3 ) e, sensor noise • an unbiased white noise signal due to sensor error. 
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Putting this in vector fonn: 

S = Strue + B + E S SN 
where, 

S is an 8x1 column vector [Sax& S Ay f, S Bxf , S Byf> Sa», S Ays , S Bxs , S Bys ] T 

Strue ^ an 8x1 column vector: [SAxCtrue? SAyfmse, Sjixfvrucj Sgyf^ue, S Axs tma^ ^>Ays true> S Bx s iruuj S B ys true] 

B is an 8x1 column vector [pi, p 2 , fh, p 4 , $5, p6> $n> Psl T 

Essn is an 8x1 column vector [sssi, S S S 2 > b S s3> e S s4, Bsss, £ss6> e S s 7 > £sss] T (“ssn” indicates Strain gauge Sensor 
Noise) 

Assume that S^, the actual force at the strain gauge, can be directly calculated from the net forces and torques on 

the rotor, [F x , F y , T*, T y ], accounting for the rotation transformation and assuming random force noise (actual forces 
at the strain gauges, such as vibrations, that do not result in a net force or torque on the rotor). Putting this m 
equation form: 

Strue " T G F nct + E S JN 
where, 

T is an 8x8 transformation matrix that accounts for the angle of the rotor relative to the fixed gauges. It is a 
calculable direct function of \\i. The upper right and lower left quadrants are all zeros. The lower right quadrant is 
the identity matrix (since the spinning geiuges do not require transformation). “T” is for “transformation." 

G is a constant 8x4 matrix that accounts for the strain gauge locations in the transformation from F net to 3^. It 
contains elements such as z Bs / (z Bs z As ). “G” is for “geometry ” 

V Mt is the 4x1 vector of net forces and torques on the rotor, [F*, F y , T* T y }. 

Esfn is an 8x1 column vector' [e SF i, e SF2> SsF 3 > e S F 4 > b S f5> £sf6> 6sf7, Bsfs] T (“sfn” indicates Strain gauge Force 
Noise) 

To first derive G, assume the rotor angle, \|/. is zero so T is the identity matrix In this case, 

[ Strue ~ GF ne t 


SAxftme 

On 

G12 

g 13 

Gh 


S Ayf true 

G 2 i 

G 2 2 

G 23 

G24 

Fy 

SBxftrue 

G 3 i 

g 3 2 

g 33 

G 3 4 

% 

S By f irue 

G41 

G42 

G43 

G44 

A 

Saxs trite 

Gsi 

G52 

g 53 

G54 


Sajts true 

Q*1 


g 63 

G64 


SBxstrue 

G71 

G72 

G73 

G74 


S;Bys true 

Gsi 

Gs2 

G83 

Qu 



Each of the elements in G will now be identified by force and torque balance equations 

Fx ^ S^xfirue "h S B xf Uue 

Ty “ S^xftoie ZAf + Sgxfhue Z B f 

solving these two equations with two unknowns. 
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SAxftnie [”ZBf / (ZAf Z B f)] F x + [1 / (z^f ' Z B f)] Ty 

SBxflnie = [?Af / ( z Af • z Bf)] F x + [* 1 / (ZAf Z B f)] T y 

Repeating fory, 

Fy — S Ayf trye + true 

%c “ *S Ayf true ZAf "F “S B yf true z Bf 

solving these two equations with two unknowns, 

$Ayf tnie ~ [~ z Bf / (ZAf Zgf)] Fy ■+ [• 1 / (ZAf Z B f)] % x 
S B yf true - [ z Af / ( z Af z Bf)] Fy + [1 ' (?Af z Bf)] T x 

With y~0, the only difference for the spinning gauges is the different z-axis gauge locations, 

Saxs true ~ [ -Z Bs / ( z As z Bs)l F x "F [1 / (Zas Zg s )] Ty 
Sbxs true ~ [ z As ! (Zas z Bs) 1 Fx "F pi 1 ( z As Z Bs )] T y 
S A ys tnie 1=5 [~ z Bs / ( z As z g s )] Fy + [■ 1 1 (Za$ Zgs)] T x 

S B ys true [z As / (z As z Bs )] F y + [\ / (z As z Bs )] T x 
So the G matrix is: 

G = “ZBf/ ( z Af Z B f) 0 0 1 / (ZAf - Z B f) 

0 -Z B f/ (ZAf Z B f) 1 / (Z^ Z B f) 0 

z Af / ( z Af z Bf) 0 0 ■ 1 / (ZAf Z B f) 

0 ZAf/ (ZAf Z B f) l/ (ZAf Zfif) 0 

“Z B s ' (z A s Z Bs ) 0 0 1 / (z A s ZBs) 

o 'Z B s t (ZAs Z Bs J 1 > (zas z 0s > 0 

z As / ( z As z Bs) ^ 0 1 P( z As z Bs) 

0 z As ' ( z As z Bs) 1 ' ( z As Z Bs ) 0 

Now allowing fory # 0, find the transformation matrix T that performs, 

Stole 5=5 F Suruel^ *0 

The spinning gauge forces are independent of the rotor angle, so the lower right quadrant is a 4x4 identity matrix 
There is no coupling between spinning and fixed gauge forces, so the upper right and lower left quadrants are all 
zeros. The upper left quadrant performs the following transformation, 

SAxftnie™ cosy -Siny 0 0 SAjcfttuelvy^Q 

S Ayf true Siny cosy 0 0 

S B xftnie 0 G COSy -Siny S B xftniely*0 

Ssyftrue 0 0 Smy COSy S B jf buelvj/ =0 

So the T matrix is, 

T « cosy -siny 0 0 0 0 0 0 

sin\|/ cosy 0 0 0 0 0 0 
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0 

0 

cosy 

-siny 

0 

0 

0 

0 

0 

0 

siny 

cosy 

0 

0 

0 

Q 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

T l will be needed at each sample, so rather than calculate it numerically at run-time it is calculated analytically here 
At each sample cosy and siny need to be calculated once only 

T 1 ~ cosy 

siny 

0 

0 

0 

0 

0 

0 

-siny 

cosy 

0 

0 

0 

0 

0 

0 

0 

0 

cosy 

siny 

0 

0 

0 

0 

0 

0 


cosy 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 0 

It is simple to check that T 4 T - 1. 

0 

0 

0 

0 

0 
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Now that G and T have been derived, equations *** and *** are combined to yield, 

S - T G F net + B + E S fn + B SS n 
R earranging, 

S B * * T G F net + (E sfn + E S sn) 

T l ( S B) = GF net + T 1 (E SFN + E SS N) 

This is now in the standard form for a least squares problem; y = Ax + e, where T' l (S-B) is y, G is A, F net is x, and T 
^Esfn + E$sh) is e (still white noise). The least-squares solution for this equation is x - (A T A)" 1 ' A T y. The former 
representation, with (S-B) as y and TG as A, could have been used, however, by pre-multiplyingbyT 1 the A matrix 
remains constant, so it (and (A T A)' 1 A T ) does not need to be re-calculated at each update. The least-squares 
minimization will be derived here for this problem. 

Problem statement 

Assume a system governed by the above equation. T 4 and G are known. 8 S is measured, but is corrupted by the bias 
and noise terms as shown. Find F llet A that can reproduce (T' i (S-B)) A by the equation (T l (S-B)) A = G F ne{ A where the 
error ((T 4 (S-B)) A (T 4 (S-B))) 2 is minimized. The idea is that by finding F net A that minimizes this cost function will be 
close to the actual F net , and this problem formulation is mathematically easy to solve. 

minimize overF n * A J - ((T 4 (S-B)) A (T 4 (S-B))) 2 


*T 4 is based on measurements of y, but this is highly accurate compared with other measurements G is based upon 

measurement of the strain gauge locations, which is assumed to be highly accurate. Additionally , the assumption that 

the system is governed by the given equation is an important one. Every effort has been made to account for all 
effects in the model, such as axis misalignment and sensor biases, but there are sure to be some effects that remain 
unaccounted for 
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J = (GF„« (T-’(S-R))) 2 
the cost! is minimized when dJ/3Fnet A = 0 

9J/9Fnet A » 2 c*/9Fnet A (G F nBt A (T 4 (S-B))) t (G F* (T 4 (S-B))) 

03/SFnet' = 2 G T (G Vj' (T l (S-B))> 
setting the derivative, dJ/dFnef ~ 0, 

0 = 2G t (GF m a (T‘(S-B)) 

0 = G t G F,,/ G T (r‘(S-B)) 

G t GF m T = G t (T 4 (S-B)) 

(G t G)' 1 (G t G) F„ a = (G t G) 4 G t (T'(S-B)) 

F„ c1 a =((G t G) 4 G t )(T- 1 (S-B)) 

Summarizing, 

F M >rT J (S-B) 

The term, ((G T G)' 1 G T ) is a fixed 4x8 matrix that is a function of the fixed system geometry only (strain gauge 
locations, etc.). It is calculated once only, and renamed V where I = ((G T G)‘ l G T ). Some of the elements of f are 
always zero, as indicated below. Also, 

r 22 = Tm t 24 ~ r i3, r 26 = r 15, r 2 $ = r n, r 32 = -r 4 i, r 34 = -r 43 , r 3 $ = -r 4 5, r 3 g ~ -r 47 . 

T is a function of ip ; and is calculated at each sample. This is an 8x8 matrix, but the inverse is performed 
analytically, so only eight terms need to be calculated at each sample. 

S is the 8x1 vector of measurements resulting from the filtering and combination of strain gauge pairs. 

0 is the 8x1 bias vector that will be calculated in the following identification, 
svp = simp. exp s cosxp in the equation below 


Fx A = 

r i, 

o r 13 o r„ 

0 r t7 

0 

C\| J 

sxp 

0 

0 

0 

0 

0 

0 

S^xf 

p. 

K 

0 

F 22 0 

r 24 0 

r 26 0 

T 2 8 

-sip c xp 

0 

0 

0 

0 

0 

0 

$Ayf 

Pz 

< 

0 

r 32 0 

r 34 g 

r 36 0 

T 3 g 

yj 

0 

0 

exp 

S\|/ 

0 

0 

0 

0 

SBxf 

P3 


Fu 

0 r 43 0 r 45 0 r 47 

0 

0 

0 

*s\p exp 

0 

0 

0 

0 

SByf 

p4 







0 

0 

0 

0 

l 

0 

0 

0 

Saxs 

p 5 







0 

0 

0 

0 

0 

1 

0 

0 

Sajb 

p* 







0 

0 

0 

0 

0 

0 

1 

0 

Sbxs 

P? 







0 

0 

0 

0 

0 

0 

0 

1 

S Bys 

Pb 


This concludes the identification of the net forces and torques on the rotor, F net A = [F X A F y A T X A T y ^\, based on the 
strain gauge measurements. 

3.5 Calculation of forces due to imbalances alone by subtracting out the effects of the counterweights. 

The net forces and torques on the rotor are due to the sum of the imbalance forces and the counterweight forces: 

Fnet 88 T nellB + TnelCW 

F^ has been estimated based' upon the strain gauge measurements, and F lie tcw has been calculated based upon the 
counterweight coordinates . Subtracting the two leaves F„ et rB> the forces and torques due to the imbalances. 
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expanding, 

FjietlB = r T* l (S-B) Fnetcw 
■where, 

F* A - Tn 0 F u 0 r 15 0 r„ 0 

Fy 0 Y%2 0 r 24 0 P 26 0 p 28 

% 0 F32 0 F34 0 p36 0 P38 

Xy F41 0 P 43 0 P45 0 1^47 0 


o\p 

s\y 

G 

0 

0 

0 

0 

0 

SaxT 

P. 

-SV| / C\\l 

0 

0 

0 

0 

0 

0 

^Ayf 

Pa 

0 L 

0 

C\|/ 

Sty 

0 

0 

0 

0 

Sfixf 

p3 

0 

0 

-Sty c\| / 

0 

0 

0 

0 

Sfiyf 

P4 

0 

0 

0 

0 

1 

0 

0 

0 

Saxs 

Ps 

0 

0 

0 

0 

0 

1 

0 

0 

S Ays 

Ps 

0 

0 

0 
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0 

0 

1 

0 

Sbxs 

p7 

0 

0 

O 

0 

0 

0 

0 

1 

SBys 

p« 


FxCW- %CW® 2 XpCW 

Fycw = m P cw<h 2 ypcw 
T x cw= -m p cwgy P cw 
t y cw = mpcwgXpcw 


+ nipcw^ypcw + 2m pC w©v >p cw m pC waxpCw 

^pcw^xpcw 2m p cw©v XI) cw nvcw^cw 

2mccwG> 2 /ycCW + 2m cC wa^ cCW + 4tnccw®/v X ccw + 2mccw^>cCW 

+ 2mc C w© 2 ^cCw + 2m cC waiy c cw + 4raccw®^ycCw 2m cC \v^ixccw 


3.6 Identification of the model-imbalance parameters based upon these net forces and torques. 

From the previous step, 

FnetIB " P T _1 (S-B) F netC W 
- From Section 3.3, 

F netIB = O 0 

Combining these equations, a least squares minimization problem can be solved to obtain the model-imbalance 
parameters, 0 * [x^, y p , x c , y c , •Vxp> Vyp, Vjx, Vyc, ^xpj a xc , z mi|>x> Z m $y] 

F liet iB will actually have a noise signal added to the forces and torques resulting from the imbalance parameters, so, 
FneflB = O Q + Enpn 

E^fn is an 8x1 column vector [bnfu Sm> Znn, £nf 4 » ^nfs, Enf6> £uf7, ^n ] T C nfn” indicates Net-Force Noise) 

where is the “Net-Force noise” signal Earlier noise vectors E S fn and Essn were for “Strain gauge Force Noise” 
(forces at the strain gauges that did not result in net forces) and “Strain gauge Sensor Noise” (sensor errors not 
resulting from actual forces). 

Combining equations *** and *** 

P T ^(S-B) F n etcw = O 0 + Enfn 

T T ' l S- P T _1 B Fnetcw * $ 0 + Enfh 

F T -1 S Fnetcw — ® G + P T "* B + Enfh 
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r T S F^cw = [®. r T 4 ] [6, B}. + Exm 

This is in the same form as the standard linear least squares minimization problem, where (f and fi A will be estimated 
The left side of the equation (the usual place for the vector of measurements) can be calculated directly, resulting in a 
4x1 vector 

The solution to this equation is: 

ie A , b a ] = (to, r t T [®, r t 4 ]) 4 [®, r t 4 ] T (r t 4 s - F netcw ) 

where (summarizing from throughout this document), 

0 A is the estimated model- imbalance parameter vector : 

6 A = [Xp A y P " tC y* % v* v». V aj' a,« A V <|> y A z^yY 

B* is the estimated strain gauge bias vector. 

B A =[pr p 2 a p 3 A p« a ps A p 6 a p 7 a p 8 i T 

O is the matrix defming the forces and torques resulting from the imbalance parameters, 

$ = 

© 2 a 0 0 0 2© 0 0 -i 0 0 0 gc\{/ gsvg 0 0 

■a od 2 0 0 -2© 0 0 0 0 l 0 0 -gsi}/ gay Q 0 

0 -g 2o if ~2 ffl 2 ^ 0 0 4 g/ 0 0 0 0 2 f 0 0 gs\(/ ~gc\j/ 

g 0 2 a?f 2c tf 0 0 0 4a/ 0 0 'if 0 0 0 g cvy g s\y 


0 0 0 2 © 

0 0 ' 2 © 0 

2 a f - 2 ©V 0 0 

2 ©Y 2 a f 0 0 


g CV| / g svj j 0 

•g si}/ g 04 / 0 


2 f 0 


'if 0 


gs\(/ -gey 
g C\f / g s\j/ 


T is a fixed 4x8 matrix that is a function of the fixed system geometry only (strain gauge locations, etc.). It is 
calculated once only Some of the elements of F are always zero, as indicated below. Also, 

F '22 = Fn, i24 =ri3, r 26 = r I5 , F 28 = Fn, F32=-r 4 i, F3 4 - -r 4 3 J r 3 6 = -'r 4 5> -tv 


r - r,, o r t3 o r l5 o r 17 o 

0 Fjj 0 Fi 4 0 F 26 O 7 Fo 8 

0 r 3 2 0 r 34 0 F36 0 r 38 

r 4l 0 r 4 3 0 r 45 o r 47 o 

r = ((g t Gy 1 g t ) 

where. 


-ZBf/ (ZAf 

?Bf) 

0 

0 


1 / (Z^ • Z B f) 

0 


'ZBf / (ZAf ZBf) 

1 ’ (ZAf 

ZBf) 

0 

z Af ■ ( z Af 

ZBf) 

0 

0 


1 / (z M ■ Z Bf ) 

0 


ZAf / (ZAf ZBf) 

l/(ZAf 

ZBf) 

0 

’ZDs / (ZAs 

Zfis) 

0 

0 


1 / (z As ZBs) 

0 


'ZBs 1 (ZAs ZBs) 

i / (ZAs 

ZBs) 

0 

z As / ( z As 

ZBs) 

0 

0 


)/(ZAs Z B$ ) 

0 


z As ■ ( Z As Zg s ) 

' / ( z As 

ZBs) 

0 
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T" 1 is an 8x8 matrix with eight terms that change at each sample all simple functions of y, 

T' 1 = cosy siny 0 0 0 0 0 0 

-siny cosy 0 0 0 0 0 0 

0 0 cosy siny 0 0 0 0 

0 0 -siny cosy 0 0 0 0 

0 0 0 0 1 0 0 0 

0 0 0 0 0 1 0 0 

0 0 0 0 0 0 1 0 

0 0 0 0 0 0 0 1 

F T 4 is a 4x8 matrix with that can be analytically calculated to save run-time computation as follows. 

r T l = Tn cosy Fiisiny Tncosy F 13 siny T# 0 Fi 7 0 

Tusiny Tn cosy -F^ smy Tn cosy 0 Fi 5 0 F n 

r 4 i siny -F '41 cosy r 43 siny r 43 cosy 0 F 45 0 F 47 

r 41 cosy F 41 siny F 43 cosy F 43 siny F 45 0 T 47 0 

S is the 8x1 vector of measurements resulting from the filtering and combination of strain gauge pairs. 

S ~ [Sax 5 ^Ayfs Sbx 5 Sgyf, Saxsj S Ays , Sbxs» S Bys ] 

SAxf = filterf (Saxxi Saxc) / 2 ) 

S Ayf ~ filterf (S Ayfi S^p)/2} 

SBxf = filterf (S sjch S.b*qV2) 

Ssyf « filterf (Sejfi S By c) ! 2 ) 

Saxs = filterf (S Axsl S AxS 2 ) / 2 ) 

S Ays “ filterf fS Ays i S Ays 2)/2) 

Sbxs ~ ’filterf (Sbxsi S Bxs 2 ) / 2 ) 

Sb>s ~ filterCfSfijsi S Bj , S 2) f 2 ) 

where S Axfl , etc. are the raw strain gauge signals, and the filter function has yet to be determined. 

Fnetcw is a 4x1 vector of the counterweight forces and torques that is directly calculated at each sample, 

Fnetcw = m p cw© 2 x P cw + m p cwGty p cw + 2m p cw©Vypcw mpcw^cw 
m pC w© 2 ypCw nipcwotxpcw 2m p cw©v xp cw nipcw a >pcw 

-m P cwgy P cw 2m cC w<fl 2 6'ccw + 2m cC wa/x c cw + 4m c cwCo/v xcC w + 2m c cw^>T:Cw 

m p cwgx P cw + 2m c owco 2 &ccw + 2m c cwa/y c cw + 4mccw^j^cw 2m c ew4*cw 

This is programmed as a recursive least squares algorithm, so the estimates for 0 and B are updated at each sample. B 
is not used for anything. The last four terms in 0 are not used either It is only the first 12 terms in that are used. The 
extra parameters identified are important since they prevent the useful parameters from being corrupted by the extra 
forces that they represent. 

So, the useful 12 parameters that are identified are, 

0 ~ [x p y p ; x c . y c v xp , Vyp , v xc v iC a xp a iP a^ a^^] 
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These parameters define the position, velocity, and acceleration of thePMI and MCI. 


3.7 SG fault detection isolation and reconfiguration (FDIR) 

If a SG is lost, F n<5t can still be identified, but not as efficiently as before. The governing equation is repeated here, 
ignoring B. 

S = T G + (Esfn + Essm) 

When a SG is lost (assuming that when it fails it reads zero plus noise - not really the case here, since we’re assuming 
each SG is in a pair reading the exact opposite force), it is as if the corresponding row in the T matrix becomes all 
zeros. Alternatively, the row can be removed from the T matrix and the S vector, giving them each 7 rows. The 
choice here is to keep 8 rows. T in this case is not invertible, so the equation cannot be multiplied byT 1 as was done 
before for computational efficiency. So for the least squares solution of F nct , A = T. G, y = S, and 

FneT ^ ((G T T r T G)- 1 G t T t > S 

The FDI approach taken here is analogous to the commonly applied bank of Kalman filters. F lust A is estimated for 9 
different conditions: the case of no failures, and then once for each situation where a strain gauge would have failed. 
The residual from each estimate is then calculated as 

residuali = 8 T, G F netil A 

where ? goes from 0 to 8, 0 indicating the case of no failures. T, is theT matrix corresponding to failure i, and F nct / is 
the net force estimated assuming failure j is present. 

3.8 A MATLAJB simulation to validate the identification scheme, including sensitivity to noise in the strain 
gauge readings 

See Appendix B for the MATLAB files, Basic results confirm that noise should not be a significant problem 1 
Newton noise in the strain gauges can be tolerated easily. Gravity imbalance forces dominate these are affected by 
the distance between the strain gauge clusters, as well as the rotation speed. Also, a higher speed will reduce the 
effect of sensor noise, as the centrifugal forces are increased by the square of the angular rate, 

The recursive least squares algorithm runs at 120 Hz on a 90 MHz Pentium in MATLAB, which includes simulation 
as well as the autobalancing calculations This will slow down if more equations are added (more strain gauges), 
digital filtering is required, or a slower processor is used. However, the code could be sped up if programmed in C 

Sample screen outputs from the graphical user interface shown in the following figures, which were taken from a 
simulation in which the centrifuge was generating 1.0 g (setpoint), the counterweight motion was not enabled 
(resulting in much higher imbalance forces), and strain gauge pair #4 (non-rotating, measuring y-axis force in the 
lower plane) has failed. 
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**” upper force 
*"”*■ lower force 
— a rotor angle 


rotor speed » 0.454 rps 
artificial gravity = 0.995 g 


Figure d: Software simulation GUI - bearing forces 

This view shows the measured bearing forces, including sensor noise ( vibrations , and imbalances. 
The blue line indicates the zero-degree angle of the rotor The red and green lines are vectors 
indicating the measured forces in the upper and lower planes. Since SG # 4 has failed, the green 
vector has very little y-component (only noise). The scale factor (16) indicates that the radius of 
the blue circle is 16 Newtons, very large because the counterweight motion has been disabled in 
this simulation. 
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upper plane imbalance, counterweights sensor FD1 




actual point-mass imbalance 
4- actual mass-couple Imbalance 
ICted point-mass Imbalance 
« ID'ed mass-couple Imbalance 
residual point-mass Unbalance 
A residual mass-couple imbalance 
Mg counterweights 




Figure 10: Software simulation GUI - imbalance tracking, Strain gauge FBIR 

The state of imbalance, imbalance-identification, and imbalance-cancellation are shown in the 
left figure for upper and lower counterweight planes. The results of the strain gauge FD1 are 
shown in the upper right, indicating correctly that SG # 4 has failed. Automatic reconfiguration 
following the fault isolation enables the identification to remain very accurate, as shown. 
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strain gauge output [N] 


Figure 11: Software simulation GUI - simulation control panel 

This control panel is used to control the simulation, star ting/ stopping the rotor spin , 
counterweights, and simulated random imbalance motion. Sensor and counterweight failures are 
also controlled. The imbalance may be “driven” manually (rather than the pseudo-random walk) 
by using the mouse to drag the blue circle joystick emulator on the right. 



Figure 12: Software simulation GUI - strain gauge outputs 

The 8 strain gauge signals are shown, plotted against the rotor angle. The vertical red line shows 
the beginning/end of the most recent revolution. If the imbalance were not moving , there was no 
vibration or sensor noise, and no failures were present, the 4 rotating gauges would read constant 
values and the 4 fixed gauges would have sinusoidal values (2 pairs 90 degrees in phase apart). 
“Stator lower y” has failed, resulting in a reading of zero + noise. 

3.9 Filters to produce \j/, co, and a from the tachometer and encoder signals 
How these signals affect the identification; 

\|/ determines the angular position of the imbalance through the calculation of siny and cosy 
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© determines the magnitude of the centrifugal forces calculated. It is especially important since these will be the 
dominant forces in the absence of gravity, and co is squared in these calculations. 

a determines the magnitude of forces due to angular acceleration. These do not appear to be significant in 
comparison to the other forces. 

If significant angular accelerations were expected, it would make sense to measure the motor current for the rotor 
motor Then an estimator could be built that used motor current, angular, and angular-rate sensors along with a 
model of the rotor drive system to estimate vj/ : co, and a. 

However, angular accelerations are expected to be small or steady (during spin-up, for example). The first attempt 
will be to use simple, independent filters to produce each signal. Analog pre- filters will be used on the encoder and 
tachometer signals to prevent aliasing. Digital filters will then be used to further condition the signal, and a difference 
of the filtered tachometer signal will be used for the angular acceleration signal. 

3.10 Filters to produce F( A *b)x> F( A+B)y , F (A , B )x> and F( A . B ) y from the strain gauge signals 

Due to the presence of significant noise, filtering the force signals will be more difficult than filtering the rotor angle, 
rate, and acceleration signals. The noise in the force signals comes from two major sources: noise due to the sensor 
itself, and noise due to extraneous vibrations in the centrifuge structure. 

It may be useful to look for correlations between the sensors and subtract out the vibration signals using adaptive 
noise cancellation. This should be an excellent research project, and may be difficult. This should be pursued only if 
more standard methods do not work 

Since these signals are used for identification, and not directly for control,- it is possible to use an acausal filter to 
achieve better performance. Taking this idea to its limit, the data could be collected for the full run (forexample, 
several seconds, resulting in thousands of data points), then it could be filtered forwards and backwards in'time (for 
example, smoothing the data using the filtfiltQ function in MATLAB). A modification of this method would be to 
take data in batches, for example, 2 seconds at a time, smooth the data, and update the ID process once every 2 
seconds. To run the ID continuously, avoiding the delay that comes with batch processing, acausal filters 1 could be 
designed that use 10 samples ahead in time to calculate the measured force It may be desirable to design these filters 
adaptively based upon the actual force data. 

A possible sequence for this would' be: 

1 ) Fix the counterweights and do not move them during this sequence. 

2) Directly design reasonably good filters, probably causal ones, to filter the force signals, 

3 ) Run the identification to determine the imbalance parameters and sensor biases 

4) Assume these to be exact and calculate the “actual forces” that would result if they were indeed correct, and 
there was no noise 

5) Adapt , the filters (adding acausal taps if they’re not already there) to make the filtered force signals most closely 
match the “actual forces” 

6) Go to step 3) and repeat until the filter weights and identification parameters converge. 

The stability of this approach is not known. It’s certainly easier to try it out than to try to prove it, During the 
iterations over steps 3) to, 6), the same set of raw data could be used, but this could lead to over-filling -- adaptation 
to account for the particular data set but failure to generalize. As a final step, allowing the counterweights to move 
may produce better results. If the effects of the counterweights are found to be . significant relative to other noise 
levels, the effects of their motion could be calculated and subtracted out as part of the force filtering system. 

During implementation of the force filtering in general, this is the procedure 

1) Choose a sample rate that is fast enough to detect the signals of interest. The point mass and couple will produoe 
two superimposed sinusoids at the rotation frequency (0.5 to 1.0 Hz). A sample rate of 20 Hz or higher should 
do well to capture the form of each of these sinusoids. A sample rate of 50 or 100 Hz would be a safe bet. 
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2) Implement analog pre-filters with break frequencies below the Nyquist frequency, equal to half the sampling 
frequency So for a 100 Hz sample rate., use 2-pole or better analog pre-filters with break frequencies at 25. Hz. 

3) Collect some real data without the counterweights moving. 

4) Design causal digital filters (for example, 2-pole Butterworth at 10 Hz), but run them using the MATLAB 
filtfiltO function (resulting in double the magnitude reduction, but with zero net phase loss the benefit of acausal 
filtering). 

5) Look at how the data is changed by the filter, and tweak the filter parameters until it looks good. If this is not 
possible, look into ways to tweak other aspects of the force measurement system. 

6) Run the identification in batch mode on the filtered data (involving several revolutions of the centrifuge). Tweak 
it until the ID looks good. Tweak the filters if necessary At this point, if the imbalance can be identified well, this 
means that the information is there, and all that remains are (possibly significant) implementation details. 

7) Run the identification with 1-2 second updates (perhaps one update per revolution) to see how quickly the ID 
parameters converge. If it is much slower than 2 seconds, it probably does not make sense to go through the 
trouble of developing the continuously- updated identification system, 

8) If convergence is quick, and development of the continuously updated ID system is warranted, the procedure 
above for acausal filter design should be followed. In addition to this design issues, communications issues may 
arise. 

Once a structure has been chosen (for example, 2-second batches or acausal filtering with 10-tap leads), it probably 
makes sense to make the filtering of the rotor angle etc., fit the same structure 

4 Calculation of the desired counterweight positions based upon the identified imbalance 

The 12 parameters defining the position, velocity, and acceleration of the PMI and MCI are updated at each sample, 
based upon the calculations summarized in Section 3.6. 

0 ~ C x p yp ? Axp Vyp a xp &yp x 0 , y c v xc v lc a XG ay,/'] 

Assuming that imbalance motions (if any) lire basically random and likely to be faster than the control bandwidth of 
the counterweights, the goal here will be to move the counterweights to counteract the imbalance positions, [x p A y p 
Xc A > yiT] If this were accomplished exactly, the velocities and accelerations would match as well 

The counterweight parameters to counteract the imbalance positions are, 

x pcw ~ (m/mpcw) x p 

ypcw^ (mp/m pCW ) y p A 

Xccw ~ (itic/niccw) Xe 

Yccw~ (Wn^cw) y c A 

The equations defining the counterweight parameters are (from Section 3.3), 


x pCW = 

(Si + 83) / 4 

YpCw~ 

(8, + 84) / 4 

x cCW “ 

(8, - 83 )/ 4 

Yccw- 

(83-84 )/ 4 


Setting m p and m c to 1 kg and combining, 
(5 1 + 5 3 )/4« (l/m pCW ) V 
(5 2 + 5 4 ) / 4 ■« (1 /m pCW ) y p 
(Si - S 3 ) / 4 = (1/nVccw) x c 
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(8 2 - 84) / 4 - (l/iTiccwly/ 

Solving for [Si,S 2) 83, 84], 

5 1 « -2 ( (x p A / m pC w) + (Xc A/ m tC w) ) 

5 2 = -2 ( (y p A / m P cw) + (y c A / m cC w) ) 

83 - ’2 ( (x p A / m pC w) OC 7 m cC w) ) 

S 4 « 2 ( (y p A / mpcw) (y c A/ niccw) ) 

These are the desired counterweight coordinates to be sent to the servo control loop 

5 Servo control loop to move the counterweights to their desired locations. 

This appears to be handled already by the motor drivers, so unless they are found to be inadequate, they will be used 
directly. 

The desired CW), etc positions will be given directly to the control card. If excessive jitter results from rapidly 
changing commands, a digital low-pass filter may be added immediately before the command is sent to the motor 
control card. 

6 Development approach / Work plan 

Unlike prior work, which involved hardware development, technology development, and software simulation, the 
proposed work will focus on technology development and software simulation/testing. 

The previous work modeled the centrifuge as having a fixed spin axis, representing a hard-mounted design, and 
matching the SSRL hardware prototype (the SSRI hardware is air-bearing supported, so two axes of translational 
motion are present, and there is also slight angular motion of the spin axis - but this any deflection is opposed by a 
stiff spring and angular motion is neglected in the existing model). This simplifies the dynamics by eliminating cross- 
axis “gyroscopic” coupling, so one of the first steps will be to extend our simulation to involve motion of the spin 
axis. 

The counterweight configuration in the present NASDA design is different from that on our hardware prototype, so 
our simulation will also be updated to account for the new configuration. The use of displacement sensors, different 
from the strain gauge sensors used in our hardware, will also be incorporated into the simulation. 

Fortunately, the architecture of our autobalancing control system will allow the sensor and counterweight changes to 
be made fairly easily, The extension of the dynamic modeling to allow spin-axis wobble and translation will be 
complex, but the control system will be updated to reflect the new equations of motion. 

The spin-axis suspension system (the VIM in the NASDA design) will be modeled to represent the VIM as closely as 
we have accurate information for (e.g., if they release detailed drawings we’ll use those, otherwise, the published 
schematics will allow us to come close). Similarly with the displacement sensor modeling - we will match the 
NASDA design as closely as possible. The aGtive. damping in the VIM ts presently designed to act like a viscous 
damper so for simplicity, we will model these dampers as passive elements (so the VIM control system does not need 
to be analyzed). 

The NASDA automatic balancing system will be modeled as accurately as possible based on published reports and 
design review materials. 

With the existing NASDA (direct autobalancing control) design and the updated SSRL (indirect autobalancing 
control) designs integrated into the simulation, we will be able to perform the following specific tests: 
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1 Experimentally evaluate the robustness of each of these methods to vibrations from the ISS. The results of 
these tests will address the first proposed research area. 

2 Monitor range of travel for the VIM under various imbalance configurations. 

3 Determine the balancing fidelity possible using both methods with the sensors present in the existing design. 
If the SSRL approach shows better results, we can quantify that, and this will address the second proposed 
research area. We can also evaluate what may be possible with different sensors. 

Sensor fault detection 

Failures in sensors, either displacement- or force-, will impact the autobalancing control, probably resulting in 
excessive vibration forces. The imbalance should still be observable even with multiple sensor failures, as long as the 
failed sensors are identified and the autobalancing system is fault tolerant. This work will leverage off other Fault 
Detection Identification and Recovery (FDIR) research by the SSRL. 

Direct vs. indirect autobalancing 

The work so far at the SSRL has used an indirect control architecture, whereas the published NASDA approach uses 
a direct architecture. The indirect approach should be able to provide better ID, and therefore, better cancellation, in 
the presence of noisy sensors. A simulation-based study would investigate and quantify tire differences between the 
two approaches. This will be more involved than die analysis proposed above, since the control systems will be tested 
more thoroughly, rather than in a few specific areas. 

Active VIM control 

The present design uses the voice coil actuators to implement viscous damping. Active control of these actuators may 
enable better vibration isolation through feed-forward control to null the spring forces (as desired) and by directly 
dealing with the gyroscopic complexity by torquing about the correct axis (rather than the present, passive approach 
that provides a torque about the wrong axis, leading to a nutation that will subsequently be damped out). 

Hardware implementation 

Although initial technology development will be developed and validated using a software simulation, implementation 
in hardware would further validate the technology The SSRL Centrifuge could be retrofitted to more closely 
approximate the ISS Centrifuge design as follows: 

• The rotating strain gauges could be replaced with soft springs/dampers to approximate the VIM 
stiffhess/damping (except with the constraint that the translational and rotational stillnesses are not 
individually selectable). This would enable the rotor to translate and wobble with respect to the bearing axis 
~ again, not a perfect approximation of the VIM, but this would introduce the gyroscopic physics. 

• With the strain gauges removed, non-rotating, non-contact eddy current displacement sensors could be used 
on a ring attached to the bottom of the rotor These sensors would emulate the VIM BDSs. 

• The counterweights could be reconfigured to more closely approximate the ISS Centrifuge design. 


Appendix A. List of variables 


Appendix B. MATLAB files 


Appendix C. Proof that model structure is general 

The imbalance' model structure has been proven to be sufficiently general to model any possible imbalance, for a rigid 
body rotating about a fixed axis. 

(hand-written notes) 

Basic approach of the proof is to state that any rigid body can be represented by an infinite collection of point masses 
Then it is shown that the resulting net forces (centrifugal, accelerative, Coriolis) summed up from these infinite 
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